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ABSTRACT 

Multiwavelength observations of Galactic black hole transients have opened a new path to 
understanding the physics of the innermost parts of the accretion flows. While the processes 
giving rise to their X-ray continuum have been studied extensively, the emission in the optical 
and infrared (OIR) energy bands was less investigated and remains poorly understood. The 
standard accretion disc, which may contribute to the flux at these wavelengths, is not capable 
of explaining a number of observables: the infrared excesses, fast OIR variability and a com- 
pUcated correlation with the X-rays. It was suggested that these energy bands are dominated 
by the jet emission, however, this scenario does not work in a number of cases. We suggest 
here an alternative, namely that most of the OIR emission is produced by the extended hot 
accretion flow. In this scenario, the OIR bands are dominated by the synchrotron radiation 
from the non-thermal electrons. An additional contribution is expected from the outer irradi- 
ated part of the accretion disc heated by the X-rays. We discuss properties of the model and 
compare them to the data. We show that the hot flow scenario is consistent with many of the 
observed spectral data, at the same time naturally explaining X-ray timing properties, fast OIR 
variability and its correlations with the X-rays, which were not possible to understand within 
the jet paradigm. 

Key words: accretion, accretion discs - black hole physics - radiation mechanisms: non- 
thermal - X-rays: binaries 



1 INTRODUCTION 

Although the black hole X-ray binaries (BHB) have been inten- 
sively studied for over four decades, many problems remain un- 
solved. Among the most debated topics are the physics of state 
transitions, the interplay between the cold accretion disc and the 
hot medium, the role of the jet, the source of rapid variability, ra- 
diative processes shaping the broadband spectrum and, specifically, 
the nature of various components contributing to its different parts. 
When addressing the latter problem, three distinct components are 
usually considered: the standard (or irradiated) cool accretion disc, 
the hot inner flow (or corona) and the jet. Their relative contribution 
depends on the spectral energy range and varies with time and can 
be assessed by performing (quasi-) simultaneous multiwavelength 
observations. 

Over the past decade, numerous multiwavelength campaigns 
have resulted in a significant progress in the field. Broadband, radio 
to X-ray, spectral energ y distributions (SEDs) for many BH sources 
were constructed (e .g. Hvnes et alj 2000l: McClintock et al. 200 ll: 
Chaty et alj l2003l : ICadolle Bel et alj l200l 1201 ll : IPurant et alj 
20091) . In addition to the spectral information, data on the fast vari- 
ability are now available in the X-rays as well as at lower ener- 
gies in the optical, infrared (OIR) and ultra-violet bands. The light 



curves at lower energ i es are significantly correlated wi th the X-rays 
jKanbach et al.l200ll:lHvnes et alj200"3ll2006„.2009b. : iDurant et alj 



l2008ll2009l : lGandhi et alj|2010l) . showing a complex shape of the 
cross-correlation function. It provides an important information on 
the interrelation between various components and gives clues to 
their physical origin. 



The radio emission in BHBs is likely dominated by the jet 
as supported by the obser ved linear polariza tion at a 1-3 per 
cent level in the hard state jCorbel et alj|200(ih and up to ten per 
cent in spatially resolved components during the transient events 
l lFenderetalJll999l : iHannikainen et alj 12000) . In addition, a rel- 
atively high luminosi ty, requiring the size exceeding the typi- 
cal binary separation i Fended l2006t). as well as th e detection of 
superluminal motion jMirabel & Rodriguezl [19941) lean towards 
this interpretation. The power-law-like radio spectrum is often at- 
tributed to synchrotron emission o f an inhomogeneous sourc e in 
analogy with the extragalactic jets jBlandford & Koniglll 19791) . In 
blaz ars, the jet is also responsible for the X-ray and 7-ray produc- 



tion dK onigl 198lVDe nner & Schhckeiseill993l : lsikora et aljl9"94l : 
IStem~& Poutanen 200^ . On the contrary, the jets in BHBs are un- 
likely to be responsible for bulk of the X-ray photons (for compre- 
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hensi ve discussion see lPoutanen & Zdziarskil2003l : lzdziarski et alj 
l2003h . 

The spectra of hard-state BHBs constitute a power-law in 
the X-ray band with a s table spectral slope and ubiquitous 
sharp cut-off at ^100 keV jGierlinski et al.lll997l : IZdziarski et"!!] 
1 19981 : llbragimov etal] |2005|) . It is broadly accepted to be 
produced by thermal Com ptonization (e.g. IPoutanenI 1 19981 ; 
IZdziarski & GierlinskilllOoi) . Additionally, a Compton reflection 
feature originating from cool opaque matter (likely the cool ac- 
cretion disc) is often detected. Its streng th is correlated with the 
X-ray slope jZdziarski et"ailll999l llOOa) . with the width of the 
iron line as well as with the quasi-periodic oscillati o n (QPO) 
frequ ency jOilfanov et al] 1 19991 : iRevnivtsev et all 1200 ll : iGilfanovl 
l2010h . These observations support a view that the X-rays are pro- 
duced in the very vicinity of the BH, in a hot flow surrounded 
by the cold disc. In this scenario, variations in the mass ac- 
cretion rate are correlated with the c ool disc truncation radius 
jPoutanen et al. I ll997l : lEsin et alj[T997l) , with the flux of soft seed 
photons that determines the spectral slope and with the reflec- 
tion amplitude that scales with the solid angle the cold disc is 
seen from the hot flow. Correlations with the QPO frequency 
are also naturally explained, if the oscillations are produced in 
the inne rmost part of the accre tion flow by Lense-Thirring pre- 
cession jingram & Pond I20T 11). Such a scenario would favour 
models, whe re seed photons for Compt onization are provided by 
the standard Ishakura & SunvaevI iwm accretion disc. However, 
the hot flow itself also produces synchrotron radiation that can 
contribute or even dominate the seed p h oton flux to the Comp- 
tonizing me d ium 1 Ghisellini et alj 1 998 ; War dzinski & Zdziarskil 
2OOOI I2OOII: IPoutaneri & Vurnj 12009.: .Malzac~Belmontl I2OO9I : 



Sobolewska et al.ll201 ll : IVeledina et alfcOl Ibl) . 

Discovery of the high-energy (MeV) tails in the hard- state 
accreting BHBs llMc Connell et al. 1994, 20021: iLing e"tai]|l997l : 
iDroulans et ai]|201(]| : Ijourdain et alj|2012l) suggests the presence 
of non-thermal particles in these systems. Such particles may be 
produced in a hot inner flow or a jet. Their association with 
the jet, however, is inconsistent with detections of even more 
rominent hi gh-energy non-th ermal tails in the soft state of BHBs 
et all 11998: Gierlinski^etal] 1 19991: Izdziarski et alj [2OOII: 
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lMcConnell"etalji2002i: jGierlihskr& Donell2003l) , when the jet is 



quenched {Fenderetal 



20041) . Thus, the hot accretion flow re- 



mains the only alternative. In the hard state, the entire X/7-ray 
hard-state spectra can be produced by hybrid (thermal plus non- 
therm al) electrons via synchrotron self-Compton (SSC ) mecha- 
nism jPoutanen & Vurnj2009l:lMalzac & Belmontl200^ . The ther- 
mal part of the particle distribution is responsible for the power- 
law-like continuum with the sharp cut-off, while the non-thermal 
particles both produce seed synchrotron photons for Comptoniza- 
tion and contribute to the MeV energies via inverse Compton pro- 
cess. Transition to the soft state can then be associated with the 
rising role of the disc as a source of seed photons, which in- 
crease Compton cooling and cause changes in t he electron distribu- 
tion f r om mostly thermal to nea r ly non-thermal jPoutanen & Coppil 
1 19981 : IPoutanen & VurmI l2009l : Iveledina et al.l l2011bl) . The SSC 
mechanism was also shown to be consistent with the peculiar 
optical variability, which in a number of BHs in low-mass X- 
ray binari es (LMXBs) is parti a lly anticorrelated with the X-ray 
emiss ion l lKanbach et"aLl I2OOII : IPurant et aO I2OO8I : ICandhi et all 
l2008h . Namely, the increasing mass accretion rate results in a 
higher X-ray and a lower synchrotron OIR emission, because of 
an increasing role of s ynchrotron self-absorption within the source 
jVeledina et ani201 lij) . 



The OIR spectra of L MXBs ofte n show an excess ab ove the 
standard accretion disc (e.g.. lHvnes et al. 20 00. 2002: G elino et alj 
In some cases the spectrum can be described by a power- 
law of index close to zero (i.e. F^, oc Such data were pre- 
vious ly explained by additi onal contribution from the irradiated 
disc i Gierlihski et al.' '2009|), dust heated by the secondary star 
jMuno & Mauerhan.,2006,) or the jet jHvnes et al.l2002l : lGallo et alj 
l2007h . We show that they also can be explained by the synchrotron 
radiation from the non-thermal particles in the hot flow. However, 
in some cases th e OIR fluxes are higher than expected from any 
candidate alone jChatv et alj|2003l: iGandhi et al.l l2010l) . suggest- 
ing contribution of at least two components simultaneously. This 
can be a reason for the complex shape of the optical/X-ray cross- 
correlation function ( I Veledina et alJuOl lal) . 

The general shape of the time-averaged X/7-ray spectrum 
of BHB can be well explained in terms of a one-zone hybrid 
Comptonization model. However, th e short-term spectral variabil- 
ity, reflected in hard X-ray time-lags jMivamoto & Kitamotol 19891 : 
lNowaketal.ll999ah and asymmetries of the cross-correlati on func- 
tion between hard and soft X-ray energy bands t Priedhor skv et al] 
1 19791 : lNolanetalj|l98ll : iMaccarone etaLll2000ir suggests that a 
number of regions simultaneously contribute to the total spectrum. 
The observed logarithmic dependence of the time-lags on photon 
ener gy can be ph enomenologically explained by spectral pivoting 
aPoutanen & Fabia n 1999_). Theoretical model capabl e of explain- 
ing th e observed timing propertie s was proposed b y Kotov et 5] 
( l200lh and further investigated in lArevalo & UttlevI l l2006h . It a~ 
sumes that the X-ray spectrum is produced in the hot corona, 
present in a range of radii, by Comptonization of the disc photons. 
The power-law slope of the locally emitted spectrum depends on 
the distance from the BH: the hardest spectra are produced in the 
innermost region due to the lack of photons from the cold disc. The 
main source of the short-term variability in BHs is believed to be 
fluctuations in the mass accretio n rate, propagating through the hot 
accretion flow l lLvubarskiill997l) . The hard time-lags thus naturally 
appear from the perturbations propagati ng from the larger d istances 
to the vicinity of the BH. The model of lKotov et alj ( I2OOII) consid- 
ers thermal Comptonization of the disc photons, but as we show 
below the results hold in the framework of hybrid Comptonization, 
with the synchrotron mechanism as the major seed photon supplier. 

It is clear, that the complete description of the timing and spec- 
tral properties of BHB requires a multi-zone model. In this paper 
we construct such a model for the extended hot accretion flow. 
This model is somewhat analogous to the inho mogeneous syn- 
chrotron models develope d for extragalactic jets l lMarscheill 19771 : 
iBlandford & Konigl|[l979l) . The difference is that in our model the 
emission originates from an inflow, not an outflow. The advantage 
of the hot flow model is that the energy input can be estimated from 
the available gravitational energy transferred to particles via some 
mechanism, while in the jet scenario the energy release is an arbi- 
trary function. 

Similarly to the previously studied one-zone models 
jPoutanen & Vumjl2009l : llVIalzac & Belmontll2009l : IVeledina et all 
l2011bh . we assume that gravitational energy is dissipated in the 
flow and is injected in the form of electrons having a power-law 
distribution. We compute the steady-state particle and photon dis- 
tributions self-consistently by solving corresponding kinetic equa- 
tions. Our aim is to understand the broadband spectral as well as 
timing properties of such a flow, and compare it to the BHB data. 

In Section [2] we give constraints on the size of the hot ac- 
cretion flow that can be derived from the observed level of the OIR 
emission. We first construct an analytical model for the hot flow as- 
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300 Rs 

Figure 1. Schematic picture of the accretion flow inner regions. Red outer component represents the multi-color cold accretion disc, truncated at 300 
Schwarzschild radii (Rs). The inner parts are occupied by a geometrically thick hot accretion flow. In our numerical model, we split the inner hot flow 
into four zones with outer radii 10i?s (violet), 30i?s (blue), lOOJJs (green) and 300iJs (yellow). 



suming power-law dependences of the main parameters on radius. 
We then proceed to the numerical model where the electron dis- 
tributions and the emitted spectra are computed self-consistently. 
In Section |3] we present the results of simulations for the model 
corresponding to the hard state of BHBs. We show that the multi- 
zone hot disc model produces flat OIR spectra resulting from syn- 
chrotron emission of non-thermal electrons at different radii. We 
model then the state transitions by decreasing the truncation radius 
of the hot flow. In Section|4] we provide a detailed analysis of the 
observational data and compare them to our model and to the jet 
scenario. We summarize our finding in Section[5] 



2 ANALYTICAL MODEL 
2.1 Geometry 

Many observational properties suggest that the standard disc in the 
hard state is truncated far away from the central object (for detailed 
descript ion and challenge s to the truncated disc scenario, see re- 
view bv lDoneetalltoOTl) . The inner part is probably occupied by 
some type of geometrically thick, optically thin, hot accretion flow, 
which is responsible for the X-ray Comptonization continuum, but 
also contributes to the longer wavelengths. One can roughly esti- 
mate the minimum size of the source that is required to produce the 
observed OIR luminosity. 

Let us first assume that OIR emission is produced by thermal 
particles. The typical temperature of the electron gas determined 
from the Comptonization cut-off is kTc ~100 keV and the typical 
IR luminosity at 1 eV can reach uLi,,i ov = lO^'^ erg s"\ Thus, 
we get the minimum size from the Rayleigh- Jeans formula 



R 



1/2 



2.3 X 10^ cm. 



(1) 



For a IOM0 BH, assumed in all calculations hereafter, this cor- 
responds to 750 Rs (here Rs — 2GM/(? is the Schwarzschild 
radius). However, many observed properties (e.g. iron line width, 
amplitude of Compton reflection, dro p of the iron li ne equivalent 
width with the Fourier frequency, see l Gilf anovll 20 1 d) suggest that 
the cold disc in the hard state is truncated at a smaller radius. Thus, 
if the hot flow lies within the truncation radius, its thermal radiation 
is unlikely to be a good candidate to produce enough OIR photons. 

However, even a weak, energetically unimportant non- 
thermal tail, in addition to the mostly Maxwellian distribu- 



tion, gives a significant rise to the synchrotron luminosity (e.g. 
IWardzihski & Zdziarskill200ll) . For example, a tail containing only 
one per cent of total particle energy increases it by a factor of 
100. Accurate calculations (see below) show that the source size 
of i? ~ (30 — 50)-Rs, an order magnitude smaller than given by 
equation ([TJ, would in principle be enough to radiate the observed 
OIR luminosity. Such a size is consistent with the above estimates 
of the truncation radius and with the typical size of the region of 
the gravitational energy release. It is worth noticing that a strong 
synchrotron emission from non-thermal particl es makes it a good 
candidate for seed photons for C omptonization jPoutanen & VumJ 
l2009l : lMalzac & Belmonlll2009l) . which implies that the SSC spec- 
trum extends from the X-rays down to the OIR band with the low- 
energy turnover determined by the maximum extent of the hot flow. 



2.2 Hard-state OIR spectra 

Spectral properties of the hot flow in the OIR band can be under- 
stood from simple analytical considerations. Let us consider the 
flow with a constant height-to-radius ratio H/R ratio, extending 
between radii i?in and i?out (see Fig. [T] for the geometry). In or- 
der to estimate the synchrotron luminosity and spectra, we assume 
that the electrons follow a power-law distribution in Lorentz factor 
ne(7) = duc/d'y = no'y^^ , starting from 7ie = 1.0. Deviations 
from the power-law at low energies do not play any role, as the 
synchrotron emission produced by these electrons is self-absorbed. 
The Thomson optical depth across the disc is assumed to follow the 
power law t{R) oc i?^". For the constant H/R this is equivalent 
to no{R) oc r/ii oc R^^^^. We further assume that the magnetic 
field depends on the distance from the BH as B{R) oc R~'^. 

Our analytical model of the hot flow is analogous to the non- 
uniform synchrotron source m odels, previously applied to the emis- 
sion f rom e xtragalactic jets ( Condon & Dressel 1973 : de Bruvnl 
197d iMarsch er 19771: iBlandford & Konigll Il979l : iKonigll 119811- 
Ghisellini et alj|T985h . For the parameters considered here, most 
of the luminosity is produced in the inner part of the source, so 
that the OIR spectrum is composed of emission components com- 
ing from different radii and is not dominated by the radiation from 
th e outer regions (these two cases are illustrated in figs 2a and 2b 
in lGhisellinietalJll985h . 

A region of the disc at a given radius emits synchrotron 
radiation, which is self-absorbed below the turn-over frequency 
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ax 




Figure 2. Contour plot of the index oqir a function of parameters /3 and 
9 for tiiree values of electron index p=2 (solid black), 3 (dotted red) and 4 
(dashed blue). 



vt- For power-law electrons, this frequency can be calcu lated as 
jRvbicki & Lightmanlll979l ; IWardziilski & Zdziarskill200lh 



P + 2 
p + 4 



Gi{p)cRrano\ 



(2) 



where vi^ = e_B/(27rmoc) is the Larmor frequency, G\ ijp) ~ 1 is 
a combination of Euler's Gamma functions (due to averaging over 
electron pitch angles), is the classical electron radius. Substitut- 
ing the constants, we get 



3 X 



((TTJloiZ) Hz, 



(3) 



where Q = 10^ Qx in cgs units. The term in brackets can also 
be written as aTnoR = f{'yt)'y^, with 7t being the Lorentz fac- 
tor of the electrons emitting at the turn-over frequency. In the later 
representation, the equation is also valid for hybrid electrons (e.g., 
Maxwellian with power-law tail), as long as the electrons emit- 
ting at the turnover frequency are in the power-law tail. The low- 
frequency cut-off for synchrotron spectrum from power-law elec- 
trons scales as 



ft oc R 



-(/3(p+2)+2e)/(p+4) 



(4) 



Again, for hybrid electron distribution one should consider scal- 
ing with radius of the power-law tail (parameter 6), which can be 
different from scaling of the total optical depth. As immediately 
follows from equation ([3](, the turn-over frequency may fall to op- 
tical and even IR wavelengths for sufficiently low magnetic field 
and/or Thomson optical depth. 

The emission at the turn-over frequency is optically thick, so 
the intensity is equal to the source function for the power-law elec- 
trons. For isotropic electrons the intensity is (averaged over pitch 
angles) 



lcG2{p) 

2^3 ' 



-1/2 5/2 



(5) 



where G2(p) — 1 (again, coming from the angle-averaging). At 
each wavelength, there is a contribution from the optically thick 
and optically thin emission from different radii. For simplicity, we 
assume that emission from each radius contributes only to its own 



turn-over frequency^ therefore the resulting spectrum of an inho- 
mogeneous synchrotron source constitutes a power-law 

i^L^ = A-K'^R^i'h ~ 2 X lO^'^i^g B^^^'^ uli^ ergs-\ (6) 

Substituting the appropriate parameter scaling and using equa- 
tion ([4](, we get the spectral index 



QOIR 



56' + /3(2p + 3) -2p-8 
I3(p + 2) +26 



(7) 



where oc z^". In a wide range of parameters /3 £ [1, 2] and 
6 G [0, 1] the resulting spectral slope lies between —0.5 and 0.5 
(see Fig.|2]l. 



3 NUMERICAL MODEL 

Analytical model developed in Section [Z2l describes only the OIR 
synchrotron spectra and is valid for purely power-law electrons. 
Such distributions may result from various acceleration mecha- 
nisms such as shock acceleration or magnetic dissipation. In the 
limit of low optical depth and weak magnetic field the electrons are 
unable to cool and the shape of the distribution stays unchanged. 
These conditions might be satisfied in quiescent state, for which the 
analytical model can be applied. During the accretion outbursts the 
matter density in the hot flow increases and the energy exchange 
and cooling processes become important, thus the initial power- 
law distribution evolves. The most important mechanisms operat- 
ing in the hot rarefied plasmas of the hot accretion flows are Comp- 
ton scattering, synchrotron emission and absorption. Coulomb col- 
lisions, bremsstrahlung, and possibly photon-photon pair produc- 
tion and annihilation. At high energies, for a continuously operat- 
ing acceleration, the steady state distribution remains a power-law- 
like, but softens because of cooling by Compton, synchrotron and 
bremsstrahlung. At lower energies. Coulomb collisions and syn- 
chrotron self-absorption efficiently thermalize particles, forming a 
Maxwellian distribution. The total particle distribution consists of a 
low-energy Maxwellian plus a high-energy tail. Such a distribution 
we call hybrid. The shape and energy content of the tail are fully 
determined by the balance between acceleration and cooling pro- 
cesses. It cannot be calculated analytically, therefore we treat this 
problem numerically. The photon spectrum emitted by the hot flow 
is computed self-consistently with the particle distributions. 

The time-scale of equilibration of electron and photon distri- 
butions for typical parameters of our model is smaller than the cor- 
responding advection time in the hot flow (see Appendix|All. Thus 
we can use an assumption that the electron and photon distributions 
are in a steady-state. We obtain them by solving the relevant kinetic 
equations. 



3.1 Model set up 

We represent the inner accretion flow with a number of separate 
regions (zones), embedded one into another (Fig. [T}. Each zone 
i has size Ri+i — Ri of the order of its height Hi ~ {Ri+i + 
Ri)/2. We assume that the inner hot flow corresponds to some type 
of rad iatively inefficient accretion flow (see review in iKato et al] 
Il998l) . In such flow, the radiative loss rate per unit area scales with 



^ Precise calculations of additional contribution from optically thin parts 
result in a slightly different normalization, while the spectral slope remains 



the same (see lMarscheJI 19771) . 
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Figure 3. Photon spectra (upper panel) and electron distributions (lower 
panel) for the hard-state model with initial electron injection index Finj = 
2.5. Electron momenta z = ^7^ — 1 are measured in units of nicC. Other 
parameters are listed in Table [T] The lines coiTespond to zone 1 (long- 
dashed), zone 2 (dotted), zone 3 (dashed) and zone 4 (dot-dashed). Sum of 
the components is shown with solid line. Red short-dashed line shows the 
slopes from analytical approximation. For further details, see Section [X2l 



radius as Qrad oc R~^^^ , and the electron density scales as p oc 
R~^^^, thus the overall scalings for the luminosity and Thomson 
optical depth of i-th zone are 



R, + i 

LiX [ Q,^i{R)RdR oc 



/Rl 



V Ri 



Ri+i 



p{R)dR oc 



fRl /RT 



(8) 



(9) 



The scaling of the magnetic field w ith ra d ius is model- 
dependent (e.g., [Sh admehri & Khajenabl l2005l : iMeiej l2005l : 
lAkizuki & Fukud 120061) . Here we assume that the magnetic pres- 
sure and the radiation pressure are equal throughout the flow, from 
which we get B cx j^- s/^ -p^^g letter scaling is the sam e as in the 
accretion flow model of IShadmehri & Khaienabil ( I2OO5I) . 

The electron heating is simulated as a power-law injection 
with slope Finj, extending between the Lorentz factors 1 and 
10'^. Net energy input in i-th zone equals to the luminosity Li. 
This energy is redistributed between the particles (electrons and 
positrons) and photons in processes of synchrotron emission and 
self-absorption, Compton scattering. Coulomb collisions, pair pro- 
duction and annihilation and bremsstrahlung emission. The domi- 
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Figure 4. Same as Fig. [3] but for the electron injection index Finj = 3.0. 
X-ray spectra are harder in this case, while the OIR spectrum is softer. 



nating cooling regime for a specific electron Lorentz factor depends 
on the parameters of the region: size, magnetic field and optical 
depth (relevant scaling is given in Appendix [Ajl. In addition to the 
internally produced radiation, we also consider soft photons from 
the cold outer accretion disc in the form of the blackbody radiation 
injected homogeneously into the system. The parameters that de- 
scribe them are the colour temperature Tcoi and the ratio of the disc 
to the i-th zone luminosity / = ^iLg^isc/ Li. Here the factor Q,i 
accounts for the fact that only a part of the disc luminosity is en- 
tering the hot flow. It is fully determined by the cold disc/hot flow 
geometry and in our case is 



Ri 

Rd,h 



(10) 



where J?d,in is the disc truncation radius, (i?i/J?d,in)^/4 is the di- 
lution factor of the i-th zone, as seen from the cold accretion disc, 
and another factor of Ri/Rd,in accounts for anisotropy of the disc 
radiation. Photon escape from each zone, giving rise to the ob- 
served emission, is modelled in terms of escape time. The kinetic 
equations for electrons and photons describing these processes are 
solved using the code described in^'unn & Poutanen (2009]). The 
total spectrum of the flow is the sum of spectra from each zone. 
The spectrum of i-th zone is computed under the approximation of 
homogeneous isotropic distributions in a sphere of radius Ri+i. 

The model has seven parameters (i) the total luminosity, (ii) 
the index of the electron power-law tail (constant throughout the 
flow), (iii) the electron Thomson optical depth and (iv) the mag- 
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Table 1. Parameters of the multi-zone hot inner flow model. 



Parameter / zone 


1 


2 


3 


4 


Ri/Rs 


3 


10 


30 


100 


Ri+i/Rs 


10 


30 


100 


300 


T 


2.5 


1.3 


0.8 


0.4 


B (10^ G) 


1 


0.25 


0.06 


0.015 


Li (lO^e ergs-i) 


6 


3 


2 


1 


kT.^f (keV) 




0.25 


0.12 


0.05 



° fcTcol is the colour temperature of radiation 
coming into the hot flow from the inner edge of the 
cold accretion disc extended up to 



netic field in the innermost regions, (v-vi) indices of their power- 
law radial dependencies 9 and and (vii) the hot flow size. 

3.2 Hard state 

In the hard state, the hot flow extends to large radii > 100J?s and 
the role of the soft photons from outer cold accretion disc is negli- 
gible. Therefore, we neglect them in the simulations and consider 
only the emission from the hot flow. We take the total luminosity 
of the flow of L = lO^^I/Edd (i-Edd IS the Eddington luminosity), 
the Thomson optical depth of the innermost zone i = 1 of r = 2.5, 
typically foun d from the X -ray/7-ray data (e.g. IZdziarski et alj 
[l998; Fronter a et alj|2001bh . and magnetic field B = 10® G. The 
radial dependencies of the parameters are given in Sect. 13.11 and 
listed for each zone in Table[T](first five rows). The results of sim- 
ulations are shown in Fig.[3]for the injection slope Finj ~ 2.5 and 
in Fig.Hfor Finj = 3.0. 

Simulations show, that larger zones generally have softer spec- 
tra. The main reason is that the outer zones are more transparent 
to the synchrotron radiation, which increases the ratio of the syn- 
chrotron to the thermal Compton luminosities. At the same time 
we see that the equilibrium electron temperature grows with radius 
from approximately 40 keV up to 180 keV. This is caused by a 
significant drop of the optical depth r in the outer zones, with a rel- 
atively slow change of the Compton j/-parameter. The X-ray spec- 
trum of the outer zones is dominated by thermal bremsstrahlung, 
because its role relative to Compton cooling <x t R/L grows lin- 
early with radius (eq. (22) in iVeledina et alj201 Ibh . The combined 
spectrum of all zones has a concave shape, exactly as observed 
jlbragimov et al.ll2005l) . 

The high-energy tail above a few 100 keV is dominated 
by Comptonization produced by the non-thermal electron tail. At 
Lorentz factor above 20, the tail has a power-law shape correspond- 
ing to index p = Finj + 1 due to synchrotron and Compton cooling. 
At intermediate 7, the distribution is curved, because of the large 
role of Coulomb collisions which pr oduce equilibrium dist ribution 
with index p = Finj - 1 (eq. (12) in IVeledina etallbOl Ibl) . 

The OIR spectrum is produced by a combination of syn- 
chrotron self-absorption peaks from different zones. The outer 
zones dominate at lower frequencies. The low-energy cut-off is de- 
termined by the size of the largest zone, which for R = 300i?s is at 
0.2 eV. At even lower energies, the spectrum is oc u^^'^. Above 
10 eV, emission from all zones is optically thin and is dominated by 
thermal Comptonization of seed non-thermal synchrotron photons. 

The X-ray spectra are harder for Finj = 3.0. This is a di- 
rect consequence of the softer equilibrium electron distribution, 
which results in a lower synchrotron luminosity and larger Comp- 



ton y-parameter. The OIR spectra are, however, softer in this case, 
because a different slope of the electron non-thermal tail results 
in a relatively low normalization of the electron distribution and, 
respectively, in the weaker synchrotron emission from the inner 
zones. 

Many of the numerical results can be understood from the an- 
alytical model if one approximates the electron distribution by a 
power-law in the energy range, where electrons emit close to the 
self-absorption frequency. For our simulations these electrons have 
Lorentz factors 7t ~ 10. 

In the case of injection slope Finj = 2.5 the optical depth of 
the power-law electrons scales the same way as the total optical 
depth, with index 6 = 1/2 (see Fig.jSjj). The slope of the elec- 
tron distribution is approximately p = 2. Putting these parameters 
(with P = 5/4) into equation (O, we get qqir = —0.13, in good 
agreement with the numerically computed slope (see Fig. [3^). 

For softer electron injection Finj = 3.0, we find that the opti- 
cal depth at the Lorentz factor 7t = 10 is nearly constant for every 
zone (see Fig.|4j3), thus for analytical approximation we take ^ = 0. 
The average electron slope at this Lorentz factor is p ~ 11/3. 
Putting these coefficients into equation Q, we get a = —0.34, 
also in good agreement with the computed spectrum (Fig.|4ji). 

The turnover of the synchrotron spectrum from each radius is 
given by equations ^ and For the case with Finj = 2.5 we 
substitute parameters of zone I: R = WRs, B = 10*^ G, r = 
10^^ (Thomson optical depth of the high-energy tail), and power- 
law slopes /3 = 5/4, 6 = 1/2 and p = 2 we obtain the scaling 

--^°^^(i(i;)"«- 

The similar scaling can be obtained for Finj = 3.0. In order to 
estimate the synchrotron luminosity we substitute the calculated 
turn-over frequency into equation l[6) 

/ON -7/8 

^^^-2xl0-(^j ergs- (12) 

which is consistent (within a factor of two) with the values obtained 
in precise numerical calculations. 

As we mentioned above, the precise values of the minimum 
and maximum Lorentz factors of the power-law electrons do not 
affect much the resulting spectra. The main model parameters (see 
Sect. 13. It can be constrained by the data. The first four parameters 
can be obtained from the X-ray luminosity and spectral slope, the 
cut-off temperature and the slope of the 7-ray tail. The final three 
parameters can then be extracted from the OIR data: the turn-over 
frequency (equation [Sj, spectral slope (equation [7} and the lumi- 
nosity (equation[6l(. 



3.3 State transitions 

A generally accepted scenario for the hard to soft state transi- 
tion involves th e motion of the cold accretion disc towards the 
compact object ( |Poutanenetal]|l997l : lEsin et alj[T997l Il998h . In 
this case, the role of the disc increases and it gradually replaces 
the synchrotron as a source of seed photons for Comptonization. 
We simulate this action by replacing the spectrum in the corre- 
sponding zone of the hot flow with a multi-c olour blackbody disc 
dShakura & Sunvaevlll973l ; lFrank et alj|2002h of an appropriate in- 
ner radius. We take the disc truncation radius i?d,in equal to the 
outer radius of the largest zone of the hot flow and we keep the 
outer disc radius at i?d,out = 3 x 10* i?s- 
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Figure 5. Upper panels: geometrical evolution at state transition. Principal components are standard accretion disc (red) and inner hot flow: zone 1 (within 
lO-Rg, violet), zone 2 (within SO-Rg, blue) and zone 3 (within 100i?s, green). Lower panels: spectral evolution at state transition. Contribution of different 
zones are marked with lines: zone 1 (long-dashed), zone 2 (dotted), zone 3 (short-dashed) and thin accretion disc (three-dot-dashed). Colour coding is the same 
as in the upper panels. The inner radius of the truncated accretion disc changes from (a) lOO/Jg through (b) SOiJs ('^) lO^S. replacing the corresponding 
zones of the hot flow. 



The additional seed photons for Comptonization are mod- 
elled by the injection of blackbody photons with temperature cor- 
responding to the colour temperature of disc inner radius: 

— (i^)lii)""('-/D"-<.. 

(see Table [Til. Given that the tra nsition occurs at almost constant 
luminosity (see iDone et alj2007h . we assume the luminosity, mag- 
netic field and Thomson optical depth of each hot flow zone remain 
the same as in the hard state (see Table[TJ. The relative contribution 
of the cold disc and the hot flow in the observed spectrum depend 
on the inclination, which we take equal to 60 degrees. The resulting 
spectra are shown in Fig.jS] 

The total spectrum is now composed of synchrotron and 
bremsstrahlung photons, as well as Comptonized synchrotron and 
disc radiation. The transition between the later two emission spec- 
tra is reflected in the overall spectral curvature at ~0. 1 keV in the 
spectra of the largest hot-flow zone (see Fig. O. In this zone the 
disc is the dominant source of seed photons. The cool disc lumi- 
nosity grows when it moves towards the black hole and its photons 
are much more energetic than those provided by the synchrotron 
mechanism, thus the X-ray spectrum softens as the transition pro- 
ceeds. At the same time, the outer regions of the hot flow collapse, 
leading to a dramatic drop in luminosity at ~0.1 eV. The (low) 
cut-off frequency increases and the OIR spectrum becomes harder 
(see Fig. [6](. Relatively small changes occur around E ^ IQ eV. 
Fig.HJb) shows the spectra in more details. Here we see that the 
pure hot-flow spectrum below the cut-off in the IR band is a power- 
law with index a = 5/2 corresponding to the optically thick non- 
thermal synchrotron. Once the truncation radius decreases to 10i?s, 
the spectrum becomes a combination of the disc and non-thermal 
synchrotron, resulting in a softer spectrum with a ~ 20 

Hence, one would expect fast change in the luminosity at OIR 
wavelengths, while (almost) no change in the UV. For instance, if 
one observes the collapse of the 100i?s zone while the 30iis zone 
is still present, there will be huge changes at ~0.5 eV, while not 

^ Further softening of the spectrum is expected if the cold disc penetrates 
into the hot flow forming a corona-like geometry. This results in additional 
cooling by the disc photons and reduction of the Compton y-parameter. 



so significant changes at ~2 eV. The opposite is expected during 
the soft-to-hard state transition: when the disc recedes, the hot flow 
occupies larger and larger radii and its synchrotron luminosity in- 
creases earlier at shorter wavelengths. 

Fig. [7] illustrates possible spectral features appearing for dif- 
ferent sizes of the outer disc radius. We see that the hot flow com- 
pletely dominates the spectrum below ~10 eV if its size is larger 
than 100J?s ■ In this case, the exact value of J?d,out does not play 
any role (unless reprocessing in the outer disc starts to be impor- 
tant). The largest changes occur for smaller truncation radius and 
large i?d,out ~ 10^i?s. For such large discs (see Fig.|7fa)), radi- 
ation in the far IR is dominated by the Rayleigh- Jeans part of the 
spectrum from the outer cold disc. The UV radiation is mostly pro- 
duced at the inner disc edge. Synchrotron from the hot flow is still 
important in the optical. 



4 COMPARISON WITH OBSERVATIONS 

In the present work we considered an inhomogeneous hot accre- 
tion flow model for the broadband spectra of the accreting black 
holes. In the hard state, when the standard cold disc is truncated at 
a large radius, the central hot region is radiating mostly via thermal 
Comptonization of the non-thermal synchrotron photons. Hot flow 
extending over a large range of radii produces a power-law-like flat 
spectrum in the OIR range. 

In the soft state, the cold disc moves in, brightens and takes 
over as a source of seed soft photons for Comptonization. This ef- 
fectively reduces the role of synchrotron radiation in electron cool- 
ing. At the same time, the reduction of the size of the emitting 
region results in an increase of the synchrotron self-absorption fre- 
quency, making the synchrotron emission in the IR band negligible. 

The scenario considered in this work is capable of reproducing 
the broadband spectra from the IR to the gamma-rays of black holes 
in all spectral states. However, the spectral data alone are not capa- 
ble of distinguishing among various models and they have to be 
considered together with other sources of information (e.g. timing 
and polarization). Below we will discuss in details various proper- 
ties of the developed hot flow model and compare them to observa- 
tions. We also compare our model to the popular jet scenario. 
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Figure 7. Same as in Fig.[6p, but for outer disc radii (a) Rd,out = 10 i?s and (b) R.d,out = 300J?s- 



4.1 Hard state 

4.1.1 X-ray spectrum and variability 

In our model the X-ray spectrum is dominated by the Comptoniza- 
tion continuum from the innermost zone where most of the grav- 
itational energy is dissipated. In the X-ray range it can gener- 
ally be described by a power-law. Outer zones of the hot flow 
have softer spectra, because of a larger role of non-thermal syn- 
chrotron and the increasing amount of the cold disc photons. The 
overall spectrum is thus slightly concave. Such spectra are consis- 
tent with those observed from the black holes. For example, the 
best studied black hol e, Cyg X-1, clearly has a concave spectrum 
( lFronteraet"aLll2001al) that can be fitted with two Comptonization 
continua l lbragimov et al.ll2005l) . 

A larger contribution from the outer zones to the soft X-rays 
should be reflected in the variability properties. Assuming that vari- 
ability is produced by pro pagation of fluctuations in the mass a ccre- 
tion rate through the disc |Lvubarskii| [l997l J Kotov et al ] |200lh , we 
expect an increase in the variability amplitude for higher photon 
energies at higher F ourier frequencies, which is indeed observed 
jNowak et al.l 1999ah . The autocorrelation function of soft X-rays 
in our model is expected to be wider than that of the hard X-rays, 
consi stent with what is measured in Cyg X-1 jMaccarone et alj 
I2OOOI) . The same effect is more obviously seen i n the Fourier- 
frequency-resolved spectra teevnivtsev et alj[l999l ; lGrifanov et alj 



l2000l) . which are softer and have larger reflection amplitude at 
low Fourier frequencies. This implies that soft X-rays are mostly 
produced in the outer zones of the hot flow, closer to the cold 
reflecting medium. The reduction of the equivalent width of the 
6.4 keV Fe line in the frequency-resolved spectra above 1 Hz sug- 
gests that the cold disc truncation radius is of the order of 10 
Schwarzschild radii jRevnivtsev et alJll999l : [Gilfanov et al.ll200(]l) . 
further supporting our scenario. Similarly, even larger inner radii 
of the cold disc were measured in the low-extinction black hole 
transient XTE J 1 1 1 8-1-480 ( Esin et al. 200 1 ; Chaty et al, 2001) . 

Another important finding is tha t the harder X-rays are de- 
layed with respect to the soft X-rays jNowak etal.lll999allbh . The 
large values of these hard time lags and their frequency-dependence 

can naturally be explained by spectral pivoting of a p ower- 
law-like spectrum jPoutanen & Fabianlll999l : |Poutanenll200lh . The 
spectral evolution can arise when the accretion rate fluctuations 
propagate towards the black hole into the zone with harder spectra 
jKotov et alj|200ll) , again consistent with our multi-zone hot flow 
model. 



4. 1.2 OIR excesses and flat spectra 

The OIR excesses above the standard disc spectr um were re- 
porte d in a number of so urces: XTE J1859 -I -226 favnes et al 



l2002h . XTE Jl 118-1-480 ( iHannikainen et alj I2OO0I : lEsin etal 
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Figure 8. Same as in Fig.|6] but with the additional contribution from the irradiated accretion disc. 
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SWIFT J1753.5-0 127 l lChiang et alj 12010) . V404 Cyg 
( iHvnes et alj l2009al) . In our model, the OIR spectrum con- 
sists of two components (Fig.O: one comes from the multi-colour 
accretion disc and another from the hot flow. The relative role 
of these components varies with the wavelength (Figs |6] and O. 
The disc spectrum is hard in the OIR band, while the non-thermal 
synchrotron from the hot flow is typically softer with qqir ~ 0. 
The second component thus produces an excess emission. 

In many cases the contribution of the non-thermal component 
is rather small compared to the disc, and it can be seen only as the 
IR excess. On the other hand, sometimes the synchrotron compo- 
nent dominates, which results in an almost pure power-law OIR 
spectrum. A good example is XTE 1111 8+480, where the spec- 
tral index aoiR = —0.15 was measured l lChatvetalj[2003h . This 
spectral index can be reproduced in our model, for example, with 
parameters Q = 0.5, p = 2.0 and /3 = 5/4 (Fig. [3] illustrates this 
case). 



4.13 Optical/X-ray cross-correlation 

In recent years, a number of simultaneous optical(IR,UV)/X- 
ray observations wit h high time resolu ti on were per f ormed 



jKanbach et ai] 1200 ll: favnes et al.l |2003| ; lOandhi et al] |2008| ; 
iHvnes et al.ll2009bl : iDurant et al.ll201 ll) , all revealing the intrinsic 
connection of the two light curves on subsecond time-scales. The 
computed cross-correlation functions have complicated shape with 
a dip in the optical light curve preceding the X-ray peak (the so- 
called precognition dip), together with an optical peak lagging the 
X-rays. 

The behaviour can be explained if the optical emission 
consists of two components: one coming from the synchrotron 
in the hot flow and another from reprocessed X-ray emission 
jVeledinaetaHbOUiJ) . Increase of the mass accretion rate causes 
an increase in the X-ray luminosity and affects the parameters of 
the hot flow, leading to a higher synchrotron self-absorption. The 
latter results in a drop of the optical emission, therefore these two 
energy bands appear anticorrelated. This is reflected in the negative 
CCF with the shape resembling that of the X-ray auto-correlation 
function. On the other hand, the reprocessed radiation is delayed 
and smeared, giving rise to a less sharp CCF peaking at positive 
lags (optical delay). The combined CCF has a complicated shape 
consistent with the data. From the point of view of the multi-zone 



consideration, with small increase of mass accretion rate the cool 
disc moves inwards and causes the collapse of the hot flow at large 
radii. Thus the suppression of the OIR emission with increasing 
X-ray radiation is also expected in this scenario. 

A different, one-peak structure of t he IR/X-ray cross - 
correlation function was found in GX 339-4 l lCasella etalj|2010h . 
suggesting the hot flow was not the dominant source of the corre- 
lated variability during their observations. However, it might still 
give significant contribution to the constant flux component, but 
less to the varying component, and therefore would not be detected 
in the timing analysis. As it can be seen from Figs [3] and |4l the 
zones giving major contribution to the IR wavelengths do not con- 
tribute much to the X-rays. Therefore, the fraction of the correlated 
variability coming from these regions is expected to be small, and 
anoth er source (likely the jet emission, as suggested in lCasella et alj 
I2OIOI) is responsible for the shape of the cross-correlation function. 

The optical correlation with t he X-rays was also detected in 
the quiescent state of V404 Cyg l lHvnes et alj 12009^ . while no 
clear radio/X-ray (nor radio/optical) correlation was found on the 
time-scales of hours, again suggesting the radio and optical emis- 
sion are not produced by the same component. 



4.1.4 Irradiation of the cold disc 

The X-ray radiation from the hot flow can be intercepted and re- 
processed in the cold disc. The irradiation strongly depends on the 
disc outer radius and the disc shape. The larger is -Rd,out, the cooler 
can be this emission. The more flared is the disc, the larger is the 
reprocessed luminosity, which typically is expected to give sig- 
nificant contribution to the OIR band, exceeding the vi scous disc 
luminosity at these energies dShakura & Sunvaev|[l973l) . Presence 
of t he irradiated dis c can be also reflected in the X-ray time-lags 
(see |Poutanenl2003 . and re ference therein) and in the optical/X-ray 
cross-correlation f unction 1 Veledina et al. 201 lah. It is also seen in 
the spectrum (e.g.. lHvnes et alJl2002j : lGierlihski et alj2009l) . 

For typical parameters of LMXBs with the disc size of 
10^^ cm, the X-ray luminosity of 10'^^ erg s^^ and 10 per cent 
reprocessing efficiency, the temperature of the outer disc is about 
1.2 eV. For an illustration, we have added the emission from the 
irradiated discs to our hot-flow spectra (see Fig. [8}. Following 
ICunninghamI ( Il976l) . we assumed the following dependence of the 
effective temperature on radius 



1.2 



R 



Rd 



-3/7 



eV 



(14) 
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and the outer disc radius of -Rd.out ~ 10^^ cm. 

As can be seen in Fig. [8] for typical parameters of the BHBs 
the spectrum of the irradiated disc peaks around 5 eV and has a 
Rayleigh- Jeans-like tail with a ~ 2 in the OIR. If the inner hot 
flow size exceeds 30-100 Rs, its synchrotron emission will domi- 
nate over the reprocessing below ~1 eV. We do not expect much of 
reprocessed emission below 0.5 eV, except for the large-period sys- 
tems (similar to V404 Cyg and GRS 1915-F105). Therefore, at suffi- 
ciently long wavelengths we expect the IR/X-ray cross-correlation 
function to have only the dip, with no peak. The overal spectral 
shape is complex, with a number of bumps corresponding to the 
standard disc (above ~ 100 eV), the irradiated disc and synchrotron 
from the hot flow (in the 1-10 eV range). A hardening of the 
spectrum observed i n GX 339^ at ~2 -3 eV jSuxton et alj|20l3 ; 
iDincer et al. 2012; Rahoui et alj|2012l) can be due to a transition 
from the hot-flow to the irradiated disc spectrum. 



4.2 State transition 

4.2. 1 Broadband spectral evolution at state transition 

During the hard-to-soft state transition the outer zones of the hot 
flow gradually collapse. This leads to a drop in luminosity at longer 
wavelengths preceding the subsequent drop at shorter wavelengths. 
Such behaviour wa s found in the UV data of the GX 339-4 during 
the 2010 outburst jYan & ^^120121) . The opposite is expected at 
the reverse transition, namely, the fast luminosity increase occurs 
earlier at shorter wavelengths. 

Recently, the entire transition of the BH transient XTE J1550- 
564 from har d to the soft state an d back was monitored in the V, I 
and H filters teussell etalllioT l|). Just before the X-ray spectral 
transition as well as after the reverse transition, significant colour 
variations occurred, as indicated by the rapid changes in the H band 
at almost constant V magnitude. 

In terms of our model, the observed colour change is related 
to the collapse/recover of a zone in the hot flow that is responsible 
for the H-band emission (see Section [331 l. As is known, the hard- 
soft and the soft- hard spectral transitio ns occur at different X-ray 
luminosities (e.g. lZdziarski et al .l2004h . This hysteresis most prob- 
ably is related to the fact that at the same luminosity the cold disc is 
further away from the central source on the rising phase of the out- 
burst, than on the decline. According to our model, the hysteresis 
has to be reflected also in the OIR spectra, namely the fast colour 
change should occur at a higher X-ray luminosity on the rising 
phase , than on the decline. This is indeed observed jRussell et alj 
l20Tlh . 

In the soft state, the accretion disc extends to the last stable 
orbit, leaving no possibility for the inner hot flow to exist. However, 
the corona should still be present as supported by the existence of 
the non-thermal tails (produced by the inverse Compton scattering 
of the disc photons). It may also produce synchrotron radiation in 
the OIR band, but likely at a lower level. 



4.2.2 Change of the X-ray radiation mechanism 

At luminosities above a few per cent of Eddington, there is a strong 
correlation between spectral i ndex and luminosity. At lower lumi- 
nosities the trend is reversed jSobolewska etaljbouh . Similarly, 
an indi cation of the reverse tr end was detected in low-luminosity 
AGNs jConstantin et alj2009l) . This was interpreted as a change of 
the source of seed photons for Comptonization from the disc pho- 
tons dominating at higher luminosities to the synchrotron at lower 



luminosities. The anti-correlation at L/L-Edd ~ lO^'^-lO^'^ can 
be reproduced within a t wo-temperature hot accretion flow model 
jNiedzwiecki et all 1 20121) . The whole spectral index - luminosity 
dependence is well explain ed by one-zone hybrid Comptonization 
model (see figs 7 and 12 in IVeledina et alj|201 Ibl) . The multi-zone 
consideration presented in the current paper follows the same pat- 
tern as the one-zone model, because the X-ray spectrum is domi- 
nated by radiation from the inner zone. 



4.3 Polarization 

The only indication of th e X-ray polarization fro m BHB goes back 
to the OSO-8 satellite jWeisskopf et alj [19771) , which measured 
3. 1±1. 7 per cent linear polarization from Cyg X-1 at2.6keV. Such 
a polarization can be produced by Compton scattering if the geom- 
etry of the X-ray emitting region is a flatten ed disc-like structure 
(H/R ~ 0.2 according to the calculations of iLightman & Shapirol 
Il976l) . The number of scatterings the X-ray photons undergo de- 
pends on the electron and the seed photon temperatures. For a 
100 keV plasma, photons double their energy in each scattering, 
hence the disc photons of a typical energy of 0.5 keV would reach 
3 keV in only 3 scatterings, while the synchrotron photons emit- 
ted at 10 eV require about 8 scatterings. Thus even if the syn- 
chrotron photons are polarized, this information is forgotten, and 
the X-ray polarization is completely determined by the geometry 
of the medium. In our model, we considered the simplest case of 
H ~ R, however, the spectral shape remains the same even for the 
flatter geometry, thus the polarization measurements are consistent 
with the hot flow model. 

The polarization degree of the hot flow radiation in the OIR 
band strongly depends on the magnetic field structure. Synchrotron 
emission from non-thermal electrons can have a large polariza- 
tion degree (up to 70 per cent) in the optically thin part of the 
spectrum. A high polarization in radio and optical reaching 30- 
50 per cent is indeed observed from extragalactic relativistic jets 
jlmpev et alj|l99ll: IWills et al.lll992l : lListeJl200ll : iMarscher et all 
I2OO2I : llkeiiri et al .1201 ih . In the optically thi ck regime, polarization 
is about 10 per cent and has a differen t sign jPacholczvk & Swiharj 
Il967l : iGinzburg & Svrovatskiilll969i) . In the hot flow model, the 
synchrotron OIR power-law spectra are produced by combination 
of self-absorption peaks from different hot flow radii, therefore, 
even for a highly ordered magnetic field the polarization unlikely 
reaches the theoretical limits. For a less ordered magnetic field, ex- 
pected in a turbulent medium, the polarization levels of synchrotron 
radiation are expected to be very low. 

Optical and UV radiation from the accretion disc may also be 
polarized up to ~11.7 per cent (parallel to the disc plane) at large 
inclinations, if the opacity is dominat ed by the electron scattering 
( IChandrasekh"^ll96d ; ISobolev|[T963h . At lower energies, absorp- 
tion in the atmospheric layers of the disc reduces the polarization 
degree i lLoskutov & Sobolev|[T98llll982h and it can drop down to 
zero in the IR. A detection of li near polarization at a few per cent 
level i n the OIR bands in BHBs JSchultz et al.l2p04l : IShahbaz et alj 



l2008l : lRussell & Fendej|2008l ; IChatv et alj|201lh is thus more con- 
sistent with being produced by synchrotron radiation with rather 
tangled magnetic field, and can originate either from the hot flow 
or from the jet. In the next section we make a detailed comparison 
between these two sources. 
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4.4 Comparison with the jet paradigm 

The simplest jet model has a conical geometry with all parame- 
ters distributed as a power-law with distance an d with electrons 
having a power-law distrib ution in Lorentz factor l lMarsche^l 19771 ; 
iBlandford & Koni gill 19791) . The magnetic field can be assumed or- 
dered or tangled, but this only influences the polarization proper- 
ties. The mathematical formulation of this model is identical to the 
analytical model for the hot flow considered in Sect.|2| Such sim- 

S ilified jet model wa s applied to the broadband spectra of BHBs 
Markoff et"al]|200ll) . In a more realistic situation, the electron 
distribution would be subject to acceleration mechanisms, as well 
as coolin g processes such as Compton, synchrotron an d adiabatic 
(as e.g. in lPe'er & CasellalZOOglJPe'er & Mark"ofil2012l) . However, 
due to a complicated physics, the energy input and the accelera- 
tion efficiency throughout the jet are generally unknown, and in 
the spectral models remain ad hoc functions. This disadvantage is 
avoided in the hot flow models, as the input energy comes from the 
liberated gravitational energy, which can be estimated analytically. 
Both models, however, suffer from a large number of parameters 
and from an unspecified acceleration mechanism. 

Apart from the different physical assumptions, the jet and the 
hot flow models predict different spectral properties. In the simple 
jet model, the synchrotron spectrum consists of the optically thin 
part with spectral slope a = — (p — l)/2 and the optically thick 
part (sum of contribution from different zones) with the same slope 
as given by equation l|7]l. The low-energy cut-off is determined by 
the jet extension and falls in the radio wavelengths. In contrast, the 
cut-off of the hot-flow spectrum is related to the truncation radius 
of the disc and likely falls in the OIR band. The break energy to 
the optically thin part of the jet is determined by the extension of 
the injection zone and typically is in the IR regime, while in the 
hot flow model it is in the UV band as determined by the size of 
the inner region. The jet optically thin synchrotron is claimed to 
extend to the X-rays, while in the hot flow model these energies are 
produced in Comptonization processes. 

Let us now assess possible contribution of the jet and the 
hot flow to various wavelengths relying on the observed spec- 
tral properties. The X-ray spectra of BHBs have sharp cutoffs 
at about 100 keV that are impossible to produce by non- thermal 
synchrotron even if th e electron distribution has an abrupt cutoff 
jZdziarski etai]|2003h . The hard spectra are also difficult to pro- 
duce by optically thin synchrotron as this requires a very hard in- 
jection. Another question is then the observed low level of the X- 
ray polarization, which is in contradiction to the theoretical expec- 
tations from the optically thin synchrotron, as well as to the lev- 
els measured in extragalactic sources. On the other hand, the X- 
ray spectral properties are well explained by the (nearly) thermal 
Comptonization in the hot flow (see Sect. l4.1.TI >. 

Atop of the X-ray power-law, the Compton reflection and the 
iron line are rather often detected features. Their amplitude and 
its correlation with the underlying spectral shape strongly argue in 
favour of the small X-ray emission region and against any beamed- 
away radiation from the jet. The analysis of the hard state X-ray 
spectra of the BHB Cy g X-1 r evealed that the spectra pivot at ener- 
gies of 10-50 keV (Zdziarski e t al. 2002). If the entire IR-to-X-ray 
continuum is produced by the same region (as proposed in the jet 
model), such pivoting would predict two orders of magnitude vari- 
ations at I eV in the hard state, which are not observed. On the 
contrary. X-ray pivoting can be easily understood in the hot flow 
model, where it can be produced by small variations of the cold disc 
radius and varying injection rate of the disc photons (see Fig.jS}. 



The observed fast X-ray variability and hard time lags can nat- 
urally be understood in the hot flow model, where the X-ray emit- 
ting region is small, and the lags are related to the viscous time- 
scale of propagating fluctuations. At the same time, the lags can 
occur from Compton scattering delays within the jet l lKvlafis et alj 
[200^, however, this model contradicts th e observed narrow ing of 
the auto-correlation function with energy jMaccarone et alj|200Q) . 

The MeV tails detected in spectra of a number of hard-state 
BHBs are likely produced by the non-thermal electrons, either by 
the optically thin synchrotron emission (as in the jet) or by the in- 
verse Compton scattering (as in the hot flow). The detailed investi- 
gation of this tail showed that it can be explained by the high-energy 
end of the jet synchrotron emission, however in this case one needs 
to assume a very hard index of the electron power-law distribu- 
tion p — 1.3 — 1.6, which is in conflict with standa rd acceleration 
models and observations (see lZdziarski et al]l20I2l . and references 
therein). At the same time, the hot flow SSC easily accounts for the 
MeV tails in the hard state as well as in the soft-state, when the jet 
is quenched. 

The OIR radi ation has often been int erpreted as a jet emis- 
sion (see review bv lRussell & Fendej|2009l) . The similarities in the 
IR and radio light-c urves in microquasars, such as GRS 1915-1-105 
( lFenderetai]|l997l) . indeed favour common source of variability. 
However, this scenario meets substantial problems in a number of 
other sources. For instance, in some systems the jet broken power- 
law model, normalized to fit t he radio fluxes, significantly under- 
predi cts the optical luminosity l lSoleri et al.l2010l : lcadolle Bel et alj 
I2OI ih . Sometimes the OIR slope is different from the radio as, 
for example, in XTE Jl 1 18-1-480 the OIR spectrum with aoiR = 
—0.15 is much so fter than the radio spectrum with qr = 0.5 
jChatv et"al]|2003l) . but much harder than expected from the opti- 
cally thin jet emission. A rather low values of the OIR polarization 
(at most a few per cent) observed in BHBs is clearly much below 
the high polarization observed in the extragalactic jets. This again 
argues against strong contribution of the optically thin jet to the 
total spectrum. The hot flow scenario can reproduce the observed 
flat OIR spectra, at the same time the slopes in the OIR and radio 
do not necessarily match, as they are produced in different regions 
(inflow and outflow). 

The two models make different predictions for changes of the 
OIR spectrum during the state transitions. In the hot-flow scenario, 
the spectrum gradually hardens at hard-to-soft transition on the 
time-scales of days to weeks, coiTesponding to the typical time- 
scales on which the cold accretion disc evolves. It softens again at 
the reverse transition. On the other hand, one would not expect sys- 
tematic changes of the jet spectral slope, as any fluctuations would 
propagate through the jet on time-scales of hours, much shorter 
than the state transition. Thus, it is expected to fade/recover with 
constant spectral shape. 

In reality, the OIR can contain contributions from a number 
of components: the hot flow, the jet, and the cold accretion disc 
(likely irradiated). It is also possible that there is a dip in the mi- 
crowave band, where the transition from the radio jet to the hot flow 
occurs. The dip can be detected in the far infrared/submillimeter 
wavelengths, which were not well covered in the past. It is thus of 
high interest to observe at these wavelengths, especially with the 
available capabilities of ALMA. 
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5 SUMMARY 

The observed OIR flat spectra and the MeV tails evidence the sig- 
nificant role of non-thermal electrons in spectral formation of ac- 
creting black hole binaries. On the other hand, the commonly de- 
tected X-ray spectral cut-offs at ~100 keV can be produced only 
by thermal particles. Whether these two populations belong to one 
component or originate from completely different places is de- 
bated. We present a model, where the entire infrared to X-ray/7-ray 
continuum is produced by one component, the inhomogeneous hot 
accretion flow, present in the vicinity of compact object. The differ- 
ence from the earlier studied hot geometrically thick optically thin 
flows is that the steady-state electron distribution in our model is 
hybrid, i.e. Maxwellian with a weak high-energy tail. 

The X-ray spectra in our model are dominated by the radiation 
of the innermost regions of the hot flow. For this reason, the model 
inherit s the advantages of the one-zone synchrotron self-C ompton 
model jPoutanen & Vurnj|2009l ; lMalzac & Belm ont" 2009*) that ex- 
plains well the X-ray spectral properties of BHB in their hard state, 
such as 

(i) stable spectra with photon index F ~ 1.6-1.9 and the cutoff 
at ~ 100 keV in the hard state, 

(ii) low level of the X-ray polarization, 

(iii) presence of the MeV tail in the hard state, 

(iv) power-law like X-ray spectra extending to a few MeV in the 
soft state, 

(v) softening of the X-ray spectrum with decreasing luminosity 
below ~10~^LEdd, 

(vi) weakness of the cold accretion disc component in the hard 
state, 

(vii) correlation between the spectral index, the reflection am- 
plitude, the width of the iron line and the frequency of the quasi- 
periodic oscillations. 

We show that the multi-zone consideration allows to understand 
many other observables in the context of the hot flow model: 

(i) hard X-ray lags with logarithmic energy dependence, 

(ii) concave X-ray spectrum, 

(iii) non-thermal OIR excesses and flat spectra, 

(iv) strong correlation between OIR and X-ray emission and a 
complicated shape of the cross-correlation function, 

(v) a complex evolution of the OIR-UV spectrum during the 
state transition. 

We present relevant analytical equations to estimate the hot 
flow parameters from the OIR data (see eqs[3]|7]and|6ll. Additional 
X-ray and 7-ray data are required to find a complete parameter set. 
However, the hot flow extent can be found from the OIR data alone 
under certain assumptions (eqsll llandll2b. 

We compare the developed model to the popular jet scenario 
and show that in a number of cases the data favour the hot flow 
interpretation. We encourage future observations in the far infrared 
and submillimeter wavelengths to provide the missing link between 
radio and infrared, which would allow us to determine the contri- 
bution of the two components to the OIR emission. 
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/ R3 

We can compare the Coulomb (the radiative cooling is faster than 
Coulomb for R < 100-Rs) and advection time-scales (again, as- 
suming = 1) 

'-^^^5x10-^1 (-^)"'\ (A6) 
r \10RsJ 

Hence, for typical parameters 7 ~1-100 and r ~0.1-I the cooling 
time-scale is shorter than the advection time. 



APPENDIX A: TYPICAL TIME-SCALES 
Al Radiative versus Coulomb time-scales 

The synchrotron cooling time can be calculated as (e.g. 



1 he synchrotron cooling 
iRvbicki & LightmaiJI 19791) 



tcooi.s = ^|T-p = —— , (Al) 

|7s| 7 + 1 \3 TTleC J 

where 7s is synchrotron cooling rate. The cooling time for Comp- 
ton scattering is similar. After a little algebra one can obtain 

7.5 X 10"^ R f L \~W R \ 
tcooi.B ~ — 7 — ^ , (A2) 

where t^b ~ 3f/B/(47rf/rad) denotes the ratio of the magnetic 
and radiation field energy densities. The typical time-scale of e-e 
Cou lomb energy exchange at the equilibrium can be estimated as 



Cou lomb energy exchange at th 
(e.g. lNavakshin &Melij|l998h 



f _ 1^^2RJ^^^ 

'-COoljCoul |. lA / \^-*) 

|7Couil 3c r In A 7(7 + 1) 

where z is the particle momentum in units of rricC, 7^,^ is the av- 
erage Lorentz factor of the electrons in equilibrium and In A is the 
Coulomb logarithm. At higher Lorentz factors the cooling is deter- 
mined by radiative processes, while for lower 7 the non-radiative 
Coulomb collisions are dominant. The relative role of radiative and 
Coulomb cooling changes with radius: 

tcool.Coul ^ 5g£^7B f J_\ ( R\ ^ ^^4^ 



tcool.s "fT V ^Edd 7 V 10i?S 

where we used 7^,^ = 1 and In A — 17. Therefore, the Coulomb 
exchange rates start dominating over the radiative cooling with an 
increasing size. 



A2 Columb versus accretion time-scales 

We assume the spectra and electron distributions in each zone are 
in equilibrium, thus the typical time-scales of equilibration should 
be much less than the dynamical time-scale at a given distance. To 
estimate the latter, we conside r the properties of an ADAF flow, de- 
rived bv lNaravan& Yil ( ll994h . Taking their radial velocity approx- 
ima tion, viscosity pa rameter 0. 1 and assuming the adiabatic index 
3/2 dOuataert & Naravaail999i) . we obtain the advection/accretion 
time 
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